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ABSTRACT 

Evolutionary models for massive stars, accounting for rotational mixing effects, do not predict any core-processed material at the 
surface of B dwarfs with low rotational velocities. Contrary to theoretical expectations, we present a detailed and fully-homogeneous, 
NLTE abundance analysis of 20 early B-type dwarfs and (sub)giants that reveals the existence of a population of nitrogen-rich 
and boron-depleted, yet intrinsically slowly-rotating objects. The low-rotation rate of several of these stars is firmly established, 
either from the occurrence of phase-locked UV wind line-profile variations, which can be ascribed to rotational modulation, or from 
theoretical modelling in the pulsating variables. The observational data presently available suggest a higher incidence of chemical 
peculiarities in stars with a (weak) detected magnetic field. This opens the possibility that magnetic phenomena are important in 
altering the photospheric abundances of early B dwarfs, even for surface field strengths at the one hundred Gauss level. However, 
further spectropolarimetric observations are needed to assess the validity of this hypothesis. 
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1. Introduction 

Fast rotation is one of the most distinctive features of massive 
stars and may give rise to large-scale fluid motions in the stellar 
interior which dredge up a significant amount of CNO-processed 
material from the stellar core to the surface (Heger & Langer 
120001 Meynet & Maeder |2000l ). On the other hand, such mix- 
ing processes can transport the boron surviving close to the stel- 
lar surface, after the star has evolved off the zero age main se- 
quence (ZAMS), to deeper layers where this fragile element will 
quickly be destroyed. Therefore, the B and CNO surface abun- 
dances of OB stars are powerful probes of rotation-related mix- 
ing phenomena throughout the entire stellar interior and at differ- 
ent evolutionary phases. As such, they can be used to constrain 
theoretical models including rotational effects, such as merid- 
ional currents or shear instabilities, and to ultimately calibrate 
some poorly-known quantities, such as the diffusion coefficients 
for the transport of the chemical elements. Of particular rele- 
vance for the present work is the recent recognition that dynamo- 
generated magnetic fields could also significantly alter the pho- 
tospheric CNO content (Maeder & Me vnet 120051 Heger et al. 
12005 ). 

The evolutionary models mentioned above do not predict any 
detectable nitrogen overabundance (i.e. above 0.1-0.2 dex) in 
B dwarfs with rotation rates below ~100 km s However, it 
has now been known for a long time (but often overlooked) that 
some apparently slowly-rotating objects actually present rela- 
tively large quantities of N-enriched material at their surface, 
with overabundances reaching up to a factor three (e.g. Kilian 
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119921 Gies & Lambert I1992I I. Such an extra amount of deep 
mixing in OB dwarfs would clearly have consequences on their 
post main-sequence evolution and other related, important issues 
(e.g. stellar populations in external galaxies). 

We have recently revealed the same peculiarity in a number 
of jS Cephei stars, which have been inferred from asteroseismic 
or line-profile variation studies to be intrinsically very slow ro- 
tators. As an illustration, the B1.5-B2 subgiant 6 Cet has a true 
rotation rate of at most 28 km s"', as deduced from the detec- 
tion of a rotationally-split multiplet in MOST data (Aerts 2006|, 
yet exhibits an unexpected N excess reaching a factor about four 
(Morel et al. 2006; hereafter M06). This is puzzling when con- 
sidering, for instance, that no contamination from CN-processed 
material has been reported in the literature for some fast-rotating 
analogs in low-metallicity environments where, in addition, the 
abundance changes should be more easily detectable (e.g. sev- 
eral B1.5 subgiants in the LMC with equatorial rotational veloc- 
ities exceeding -130 km s Korn et al. I2005K To shed more 
light on this issue, here we extend the M06 study and present an 
NLTE abundance analysis of a sample of 20 dwarfs/(sub)giants 
comprising stars not established as pulsating variables, with spe- 
cial emphasis on the possible role played by magnetic fields to 
shape the CNO surface abundance patterns of early B-type stars 
during the early phases of their evolution. 

2. Target selection and observational material 

We have complemented the sample analyzed by M06 with six 
early B-type dwarfs defined as boron-poor from the analysis of 
the B HI /12066 UV resonance fine ( Brooks et al. '2002', Mendel 
et al. 2006; Proffitt & Ouigley lSOOTl Venn et al. 2002). All stars 
have boron abundances at least one order of magnitude lower 
than the meteoritic value. The studied sample includes all boron- 
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depleted B stars known to date. For comparison purposes, we 
also analyzed the boron-normal star /? Cru (ProfRtt & Quigley 
l200TT l. 

In the case of l CMa and e CMa we made use of a mean 
FEROS spectrum created by averaging 6 and 4 individual ex- 
posures closely taken in time, respectively. These data were ob- 
tained in April 2005 with the ESO/MPI 2.2m Telescope and were 
kindly provided by P. Dufton. Echelle spectra of y Ori (mean of 
two exposures) and HD 36591 were retrieved from the ELODIE 
archivefl Snapshot spectra of Ori and 15 CMa (4 and 3 expo- 
sures obtained over a timespan of about one week, respectively) 
were obtained in March 2007 with the CORALIE spectrograph 
attached to the ESO/Euler 1 .2m Telescope at La Silla. For the p 
Cephei star (3 Cru, we used a mean CORALIE spectrum created 
by co-adding a set of 49 individual exposures obtained over the 
period 21-29 March 2007. 

In addition, we discuss here previous results obtained fol- 
lowing exactly the same procedures for two slowly pulsating B 
stars Cas and HD 85953; Briquet & Morel EoOTIi, the mag- 
netic star T Sco (Hubrig et al., in prep.) and the j3 Cephei star 6 
Oph A (Briquet et al. |2007| ). It is important to emphasize that the 
performed abundance analysis is identical in all respects for all 
the stars discussed in the present paper. The abundance data pre- 
sented here thus constitute a fully homogeneous database suit- 
able for a comparative study. This point is essential considering 
the mild nitrogen excesses detected, as will be shown below. 

The majority of the stars in our sample have been the sub- 
ject of sensitive spectropolarimetric observations, which led to 
the detection of a weak-longitudinal, magnetic field in y Peg 
(Butkovskaya & Plachinda 2007), ( Cas (Neiner et al. |2003aD , 
^1 CMa (Hubrig et al. 2006 ), HD 85953 (Hubrig et al. l2006b 
T Sco (Donati et al. 2006b), V2052 Oph (Neiner et al. 2003bll 
and j6 Cep (Henrichs et al. i2007) . Very recent observations with 
FORSIA'LT not only confirm the field detection in ^' CMa, but 
also reveal a weaker field at the one hundred Gauss level in 5 
Cet, 15 CMa and e CMa (Hubrig et al., in prep.). On the other 
hand, null detections have been reported from observations with 
similar sensitivity limits for v Eri, (3 CMa, fi Cru, V836 Cen and 
e Oph A (Schnerr et al. |2006at Hubrig et al. 120061 Hubrig et 
al., in prep.). The low-sensitivity circular polarization observa- 
tions of 12 Lac by Rudy & Kemp d 19781 1 rule out the existence 
of a magnetic field at the kG level. Table [1] summarizes the ba- 
sic magnetic properties of the program stars and provides, for 
those with detections, the rms longitudinal field, {Bi), and the 
reduced of the measurements. The former is a measure of the 
intrinsic strength of the field, while the latter can be used to as- 
sess the reliability of the detection (see, e.g., Hubrig et al. 2006). 
These quantities are defined from individual measurements 
of a mean longitudinal magnetic field, (B/),, with an associated 
uncertainty cr,-, as: 



1/2 



(1) 



and 



(2) 



It should be noted that the amount of spectropolarimetric 
data collected — and therefore the robustness of the magnetic 



Table 1. Longitudinal magnetic field properties of the program 
stars. Note that the rms longitudinal field values, {Bi), are only 
indicative in case of very few measurements (e.g. 15 CMa and 
e CMa). A field is considered undetected if it does not have a 
significance of at least 3cr. 
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3.8 


1 
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38 118 


2.0 


7 




U l_-CL 


85 9 


9.5 






V En 


Not detected 








tt'^ Ori 
71 Wll 


Not observed 








y Ori 


Not observed 










Not observed 






44743 


yS CMa 


Not detected 




3,4 


46328 


^' CMa 


287 11 


478 


3,4 


50707 


15 CMa 


144 2 


19 


3 


51309 


iCMa 


Not observed 






52089 


6 CMa 


149 2 


14 


3 


85953 




102 5 


6.9 


3,4 


111 123 


ySCru 


Not detected 




4 


129 929 


V836 Cen 


Not detected 




3,4 


149 438 


T Sco 


42 35 


93 


6 


157 056 


eOph A 


Not detected 




4 


163 472 


V2052 Oph 


262 112 


1.3 


7 


205 021 


ySCep 


80 124 


8.4 


8 


214 993 


12 Lac 


Not detected 
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al. ( I2003at : [3] Hubrig et al., in prep.; [4] Hubrig et al. ( [25561 : [5] 
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field (non)detection — varies dramatically from one star to an- 
other: while few targets have been the subject of intensive, multi- 
epoch observing campaigns, others only have a handful of snap- 
shot observations available. This important point will be dis- 
cussed further in Sect.li 



' See |http : //atlas . obs-hp . fr/elodie/ 1 



3. Methods of analysis and sample characteristics 

The non-local thermodynamic equilibrium (NLTE) abundances 
of He, C, N, O, Mg, Al, Si, S and Fe were calculated using the 
latest versions of the line-formation codes DETAIL/SURFACE 
and plane-parallel, fully line-blanketed Kurucz atmospheric 
models (Kurucz4993i A standard, iterative scheme is first used 
to derive the atmospheric parameters purely on spectroscopic 
grounds: T^s is determined from the Siii/m/iv ionization bal- 
ance, log g from fitting the coUisionally -broadened wings of the 
Balmer lines, and the microturbulent velocity, ^, from requiring 
the abundances yielded by the O n features to be independent of 
the line strength. The ionic species C iii, N iii and S iii are poorly 
treated in our model atoms and only a few spectral features are 
measurable, but there is no evidence in the BO-Bl stars for sys- 
tematic discrepancies with the abundances yielded by the much 
more numerous lines of lower ionization stages (namely C n, N ii 
and S ii). All stars are slow rotators (see below) and classical 
curve-of-growth techniques were used to determine the abun- 
dances using the equivalent widths (EWs) of a set of unblended 
lines. The abundance uncertainties take into account both the 
line-to-line scatter and the errors arising from the uncertainties 
in the atmospheric parameters. Finally, we uniformly derived the 
projected rotational velocities, 0/?sin/, by comparing the pro- 
files of a set of isolated O ii lines with a grid of rotationally- 
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broadened synthetic spectra. The reader is referred to M06 for 
complete details. 



To correct for the contamination of the spectrum by the sec- 
ondary in the single line binary fi Cru, we based our analysis on 
synthetic, composite spectra assuming the following parameters 
for the companion: reff=22000 K, log ^=4.0 dex [cgs] (Aerts 
et al. 119981 ) and a microturbulent velocity, ^=5 km s \ typical 
of B-type dwarfs. We considered that the flux ratio between the 
two components is 8% in the optical band (Hanbury Brown et 
al. [1974, Popper 1968) and that the secondary has a chemical 
composition typical of OB dwarfs in the solar vicinity (Daflon 
& Cunha |2004t the value for iron was taken from M06). To 
examine the sensitivity of our results to these various assump- 
tions, we have repeated the abundance analysis after varying the 
adopted effective temperature, surface gravity, chemical com- 
position, and luminosity of the secondary within the range of 
plausible values. Namely, we assumed in turn: reff=20000 K, 
log §=3.7 dex [cgs], the abundances of all the metals enhanced 
by 0.3 dex relative to solar and a flux ratio in the optical be- 
tween the two components amounting to up to 12%, while keep- 
ing the other parameters unchanged. As can be seen in Table |2l 
taking the companion into account (regardless of its characteris- 
tics) has a negligible impact on the abundance results. The same 
conclusion holds in the case of the single line binary 9 Oph A 
(see Briquet et al. 120071) . The stars y Peg and p Cep are asso- 
ciated with faint secondaries (Roberts et al. 120071 Schnerr et al. 
I2006bl ). and their contribution was neglected. The giant tt'* Ori is 
also part of a massive binary (as confirmed by the radial velocity 
changes observed in our own dataset) with an orbital period of 
about 9.5 d (Pourbaix et al. 120041) . but the lack of information 
about the secondary forced us to treat it as a single star (Gies & 
Lambert IT992I classify n'^ Ori as B2 III + B2 IV, but no details 
are given). 

The adopted atmospheric parameters are provided in Table[3] 
The same parameters, within the errors, are also obtained for 
all studied stars when analyzing our data with the unified code 
FASTWIND (Puis et al. 12005^ and (semi-)automatic line-profile 
fitting techniques (Lefever 2007 ). Our sample is composed pri- 
marily of apparently slow rotators with low tlR sin / values not 
exceeding ~65 km s"' (these figures are upper limits in the 
case of the Cephei variables, as the lines are often signifi- 
cantly broadened by pulsations). Our projected rotational veloc- 
ities have been derived from profile fitting, but Fourier trans- 
form methods can also be used and have the advantage of dis- 
entangling different line-broadening mechanisms (Simon-Diaz 
& Herrero 120071 ). Lefever ( I2007I I used these techniques for 15 
stars in our sample and found systematically lower OR sin / val- 
ues. Although the differences are small and typically amount to 
only a few km s"', as expected considering that dwarfs suffer 
little macroturbulent broadening (Ebbets ,1979j . this bolsters our 
conclusion regarding the low rotation rate of our targets. Most 
importantly, asteroseismic and line-profile variation studies in- 
dicate that 10 stars have true rotational velocities well below 100 
km s"' or, in some cases, as low as a few km s"' (v Eri, V836 
Cen and r Sco). As detailed in Table[3] this quantity is derived ei- 
ther from: a) identifying the recurrence timescale of the changes 
affecting the UV wind lines, with the rotational period in the 4 
objects displaying phase-locked variations, orb) identification of 
the components of a rotationally split mode or theoretical mod- 
elling of the line-profile variations in the 6 pulsating variables. 



4. Reliability of the abundance results 

The CNO abundances for the complete sample are provided in 
Table m along with the NLTE literature values for boron. Most 
B-depleted stars only have an upper limit to their boron abun- 
dance and we only considered in that case the most stringent 
estimate in the literature. Table |5] presents the He, Mg, Al, Si, 
S and Fe abundances for the stars analyzed in the present study. 
Before discussing the chemical properties of our targets, we first 
assess the robustness of our results against the adopted atmo- 
spheric parameters. 

Of particular concern is the uncertainty surrounding the de- 
termination of the microturbulent velocity in the giants. As al- 
ready largely discussed in the literature (e.g. Daflon et al. 20041 
Kilianll992t Trundle et al. 120041 ). larger microturbulence values 
are usually obtained for the more evolved objects (^ and \ogg 
are tightly antic orrelated in our dataset). Furthermore, discrepant 
values are obtained depending on the chemical species chosen 
(e.g. the On features yielding systematically larger values than 
the N n lines). Fortunately, the prime indicators for an N enrich- 
ment, namely the ratios of the C, N and O abundances ([N/C] 
and [N/O]), are largely insensitive to this problem. In the case of 
e CMa, for instance, these ratios differ by a negligible amount of 
at most 0.04 dex when adopting the microturbulence estimated 
from the oxygen (^=16 km s"') or nitrogen lines (^=12 km s"'). 
The microturbulent velocity is also poorly constrained in t CMa 
because of the lack of strong ionic features in this temperature 
range, but once again this does not affect our conclusions. 

Figure[T]displays the variations of the CNO abundances and 
their ratios, as a function of the atmospheric parameters and pro- 
jected rotational velocity. The oxygen abundance becomes in- 
creasingly difficult to measure with confidence as one progresses 
towards later spectral types because of its strong sensitivity to 
errors on Teff and the intrinsic weakness of the O ii features, but 
there is evidence for lower abundance values as Teff decreases. 
The lack of similar trends, either for carbon and nitrogen, or for 
helium and the other dominant metals (except Al; see Fig|2]) sug- 
gests that this problem is most likely intrinsic to oxygen (perhaps 
related to the NLTE corrections) and not attributable to an incor- 
rect temperature scale. Because of this temperature-dependence 
affecting oxygen and the fact that carbon is a more sensitive 
indicator of core-processed material brought up to the surface 
(the predicted depletion factors are higher than for oxygen), we 
will primarily focus on the [N/C] ratio in the following (further- 
more, the uncertainties affecting the C abundances are a factor 
~2 lower; Table |4]i. 

There is also a discernible correlation between the C abun- 
dances and log g (Fig[T]). If we use the surface gravity as a proxy 
for the evolutionary status, this would interestingly point to a 
progressive depletion of carbon at advanced, stellar ages as a 
consequence of CNO-cycle burning. However, such a C deple- 
tion in the most evolved objects should be systematically coun- 
terbalanced by an N excess, which is not observed. In addition, 
the C depletion factors (based on the N excesses discussed be- 
low) are expected to be of the order of the uncertainties (~0.1 
dex) and, as such, barely detectable (there is indeed no anticor- 
relation between the C and N abundances). This, along with the 
fact that the [N/C] and [N/O] ratios show no trends with log g, 
argues against the existence of evolutionary effects in our data 
(we will return to this issue in Sect. 16. It . 

We note that uncertain, but generally low abundances are ob- 
tained for He and most metals in the case of the B2.5 II star t 
CMa (Figures [U and |2]i. These low values might indicate a gen- 
eral problem with our analysis of this star (neglect of wind ef- 
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Table 2. Sensitivity of the derived abundances of /? Cru on the assumed properties of the secondary. We quote the abundance 
differences with the values provided in Tables |4] and |5] 



Ionic species used Without companion With companion 

Areff=-2000K Alogg=-0.3 dex A log f=+0.3 dex flux ratio=12% 



He/H 


He I 


+0.005 


+0.000 


+0.002 


+0.000 


-0.003 


Aloge(C) 


Cii(+C III .14187) 


+0.03 


+0.02 


-0.01 


-0.05 


-0.02 


A log e(N) 


N II (+N 111/14634) 


+0.00 


+0.02 


-0.01 


-0.02 


+0.00 


Aloge(O) 


On 


-0.04 


+0.01 


-0.01 


-0.01 


+0.02 


A log e(Mg) 


Mgn 


+0.03 


-0.02 


+0.00 


-0.02 


-0.02 


A log e(Al) 


Aim 


+0.01 


+0.02 


+0.00 


-0.03 


+0.00 


Aloge(Si) 


Si iii/iv 


-0.04 


+0.01 


+0.00 


-0.01 


+0.02 


Aloge(S) 


Sm 


-0.03 


+0.01 


-0.01 


-0.01 


+0.01 


A log e(Fe) 


Fern 


+0.00 


+0.01 


-0.01 


-0.03 


+0.00 



Table 3. Adopted atmospheric parameters and (projected) rotational velocities. The MK spectral types are based on our T^ff estimates 
and the calibrations for dwarfs of Crowther ( 119981 ). The calibration for supergiants was used for t CMa. The 1-cr uncertainty on Teff 
is 1000 K for all stars. The quoted uncertainty on fIR sin / is the line-to-line scatter. 



HD 


Usual 


Spectral type 


Remarks" 


Teff 


logs 




QR sin /" 


Ref. 




Ref. 


number 


name 






[K] 


[cm s"^] 


[km s-'] 


[km s-'] 




[km s-'] 




886 


7 Peg 


B1.5-B2 IV 


SBl,/?Cephei 


22 500 


3.75+0.15 


l-I 


10±1 


1 






3360 


f Cas 


B2IV 


SPB 


22000 


3.70+0.15 


1±1 


19±1 


2 


55±28 (UV) 


6 


16582 


5 Get 


B1.5-B2IV 


yS Cephei 


23 000 


3.80+0.15 


i-\ 


14±1 


1 


14 or 28 (S) 


7 


29248 


V Eri 


B1.5-B2IV 


P Cephei 


23 500 


3.75+0.15 


10±4 


36±3 


1 


6(S) 


8 


30836 


Ori 


B2III 


SB I 


21500 


3.35±0.15 


8±4 


43±3 


3 






35468 


7 Ori 


B2 IV-III 




22000 


3.50+0.20 


13±5 


51±4 


3 






36591 




Bl V 




27 000 


4.00+0.20 


3+2 


16±2 


3 






44743 


^CMa 


B1.5 IV-III 


yS Cephei 


24000 


3.50+0.15 


14±3 


23±2 


1 


31 ±5 (LPV) 


9 


46328 


^' CMa 


B0.5-B1 IV 


P Cephei 


27 500 


3.75±0.15 


6±2 


10±2 


1 






50707 


15 CMa 


Bl IV-III 


P Cephei 


26000 


3.60+0.15 


7±3 


45±3 


3 






51309 


tCMa 


B2.5 II 




17 500 


2.75+0.15 


15±5 


32±3 


3 






52089 


eCMa 


B1.5-B2III 




23 000 


3.30+0.15 


16±4 


28±2 


3 






85953 




B2IV 


SPB 


21000 


3.80+0.15 


1±1 


29±2 


2 






111 123 


/JCru 


B0.5-B1 IV 


SB 1,^6 Cephei 


27 500 


3.65+0.20 


11±3 


48±3 


3 






129 929 


V836 Cen 


B1.5 V 


P Cephei 


24 500 


3.95+0.20 


6±3 


15±2 


1 


2(S) 


10 


149438 


T ScO 


09.5 V 




31500 


4.05+0.15 


2±2 


8±2 


4 


6(UV) 


11 


157 056 


eOph A 


B1.5 V 


SB l,/3 Cephei 


25 000 


4.10+0.15 


4+2 


31±2 


5 


29 (S) 


5 


163 472 


V2052 Oph 


B1.5-B2 V 


p Cephei 


23 000 


4.00+0.20 


t-i 


62±3 


1 


56 (UV) 


12 


205 021 


/JCep 


Bl IV 


SB 1, ,6 Cephei 


26000 


3.70+0.15 


6±3 


29±2 


1 


26 (UV) 


13 


214 993 


12 Lac 


B1.5 IV 


P Cephei 


24 500 


3.65+0.15 


10±4 


42±4 


1 


45 (LPV) 


14 



Key to r eferences: [1] M 06; [2 ] Briquet & Morel ( [2007]l; [3] this paper; [4] Hubrig et al., in prep.; [5] Brique t et al. ( I2007t ; [6] Neiner et al. 
I l2003ab ; [7] Aerts et al. (2006); [8] Pamyatnykh et al. ( 12504] ): [9] Mazum dar et al. ( 12006^ : [10] Dupret et al. i2004l ); [11] Donati et al. ( I2006bb : 

[12] Neiner et al. {2003b ); [13] Henrichs et al. ( 2007|; [14] Aerts ( fT996l ). 

" We only classify as p Cephei stars the confirmed candidates in the catalogue of Stankov & Handler (2005 ). 

* Note that the spectral lines are significantly broadened by pulsations in the p Cephei and SPB stars (see M06). This explains why the quoted 
QR sin i values, which reflect the total amount of line broadening arising from pulsations and rotation, are occasionally larger than the true 
rotational velocities, Q.R. 

UV: calculated from the stellar radius derived from evolutionary tracks and assuming that the rotational period can be identified with the 
recurrence timescale of the changes affecting the UV wind lines, S: from asteroseismic studies, LPV: from modelling of the line-profile variations 
arising from pulsations. 



fects being a possibility) and turn out to contribute for a large 
part to the correlations discussed above. 

Due to our inability, as others, to satisfactorily model some 
strong C n lines (most notably C ii /14267), we based our abun- 
dance determination for carbon on a number of weak lines (see 
1VI06 for the complete line list). Nieva & Przybilla ( 120071 ) re- 
cently developed a carbon model atom that allowed them to over- 
come this problem and to obtain consistent results for all lines 
they considered. A straight comparison with their results for the 
two stars in common shows that our values are lower by 0. 14 and 
0.11 dex for HD 3659 1 and t Sco, respectively. However, the net 
effect of the model atom can be better appraised by recomputing 



the abundances using their adopted atmospheric parameters and 
measured EWs. For the lines in common, we now find mean val- 
ues lower by 0.06 and higher by 0.02 dex for HD 36591 and t 
Sco, respectively. Our ability to closely reproduce the results of 
Nieva & Przybilla (2007) validates the use of our model atom. 
The systematically higher abundances they found cannot be ex- 
plained by the slight differences in the atmospheric parameters, 
but rather by the fact that their observed line profiles are ap- 
parently stronger (their EWs are on average ~15% higher than 
ours), perhaps as a result of the different continuum rectification. 

The boron abundances collected from the literature have 
been derived from the spectral synthesis of the B m /12065.8 res- 
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Table 4. Mean NLTE B and CNO abundances (on the scale in which log e[H] = 12). The number of used lines is given in brackets. We 
define [N/C] and [N/O] as log[e(N)/e(C)] and log[e(N)/e(0)], respectively. The last column indicates the chemical group (Groups 
I, II or III) to which the star belongs to (see Sect.|5]and FigO. 



HD number 


Usual name 


loge(B) 


Ref. 


log 6(C) 


log6(N) 


loge(O) 


[N/C] 


[N/O] 


Ref. 


Group 


886 


7 Peg 


2.23+0.20 


1 


8.20+0.05 (8) 


7.58+0.11 (23) 


8.43+0.28 (21) 


-0.62+0.12 


-0.85+0.30 


5 


I 


3360 


<f Cas 


<1.21 


1 


8.16+0.08 (6) 


7.97+0.13 (20) 


8.38+0.30(18) 


-0.19+0.15 


-0.41+0.32 


6 


III 


16582 


(5 Get 


1.16±0.15 


2 


8.09+0.08 (6) 


8.05+0.11 (26) 


8.45+0.26 (19) 


-0.04+0.14 


-0.40+0.29 


5 


III 


29248 


y Eri 


2.45+0.40 


1 


8.24+0.12 (10) 


7.87+0.09 (18) 


8.51+0.24(19) 


-0.37+0.15 


-0.64+0.26 


5 


I 


30836 


Ori 


<1.0 


3 


8.19+0.09 (5) 


7.54+0.15 (13) 


8.40+0.29 (24) 


-0.65+0.18 


-0.86+0.33 


7 


II 


35468 


7 Ori 


<1.01 


1 


8.11+0.09 (6) 


7.90+0.16(14) 


8.16+0.27 (12) 


-0.20+0.19 


-0.26+0.32 


7 


III 


36591 




<1.27 


4 


8.19+0.05 (8) 


7.66+0.12 (22) 


8.62+0.12 (27) 


-0.53+0.13 


-0.95+0.17 


7 


II 


44743 


/3CMa 


2.76+0.20 


1 


8.16+0.11 (9) 


7.59+0.14 (26) 


8.62+0.18 (30) 


-0.57+0.18 


-1.03+0.23 


5 


I 


46328 


^' CMa 


<1.36 


1 


8.18+0.12(9) 


8.00+0.17 (34) 


8.59+0.17 (34) 


-0.18+0.21 


-0.59+0.24 


5 


III 


50707 


15 CMa 


<1.5 


2 


8.18+0.10(5) 


8.03+0.15 (23) 


8.64+0.18 (33) 


-0.15+0.19 


-0.61+0.24 


7 


III 


51309 


( CMa 


<0.84 


1 


7.89+0.19 (4) 


7.86+0.34(15) 


8.25+0.44(11) 


-0.03+0.39 


-0.39+0.56 


7 


III 


52089 


6 CMa 


<1.62 


1 


8.09+0.12 (7) 


7.93+0.14 (26) 


8.44+0.18 (29) 


-0.16+0.19 


-0.51+0.23 


7 


III 


85953 








8.16+0.14 (9) 


7.66+0.20(11) 


8.41+0.32 (16) 


-0.50+0.24 


-0.75+0.37 


6 




111 123 


ySCru 


2.48+0.20 


1 


8.04+0.10(7) 


7.61+0.17 (19) 


8.59+0.16 (28) 


-0.43+0.20 


-0.98+0.24 


7 


I 


129 929 


V836 Cen 






8.34+0.13 (8) 


7.73+0.10(21) 


8.49+0.24 (25) 


-0.61+0.17 


-0.76+0.26 


5 




149438 


T ScO 






8.19+0.14(15) 


8.15+0.20 (35) 


8.62+0.20 (42) 


-0.04+0.25 


-0.47+0.29 


8 




157 056 


eOph A 






8.32+0.09 (7) 


7.78+0.10(23) 


8.58+0.26 (27) 


-0.54+0.14 


-0.80+0.28 


9 




163 472 


V2052 Oph 






8.21+0.07 (4) 


7.99+0.17 (10) 


8.39+0.30(14) 


-0.22+0.19 


-0.40+0.35 


5 




205 021 


jSCep 


<0.90 


2 


8.02+0.10(11) 


7.91+0.13 (19) 


8.47+0.14 (30) 


-0.11+0.17 


-0.56+0.19 


5 


III 


214 993 


12 Lac 


2.1+0.2 


3 


8.22+0.12(6) 


7.64+0.18 (17) 


8.42+0.23 (22) 


-0.58+0.22 


-0.78+0.29 


5 


I 



Key to refere nces: [1] Proffitt & Quigley ('lOOTl; [2] Venn et al. ( l2002t : [3] Mendel et al. ( l2006l l: [4] Brooks et al. f2002t : [5] M06; [6] Briquet & 
Morel J2007t ; [7] this paper; [8] Hubrig et al., in prep.; [9] Briquet et al. i2007j . 



Table 5. Mean NLTE He, Mg, Al, Si, S and Fe abundances (on the scale in which log e[H] = 12) for the stars analyzed in the present 
study. The data for the remaining objects in the sample can be found in the references in TablelH The number of used lines is given 
in brackets. For comparison purposes, the first line gives the typical values found for B dwarfs in the solar neighbourhood (M06; 
Daflon & Cunha l2004[ Gummersbach et al. |T998l Kihan-Montenbruck et al. |l994l l. 



HD number 


Usual name 


He/H 


log6(Mg) 


loge(Al) 


loge(Si) 


log6(S) 


log e(Fe) 




B dwarfs 


-0.085 




-7.4 


-6.1 


-7.2 


-7.2 


-7.3 


30836 


Ori 


0.074+0.034 (9) 


7.33: 


1:0.34(1) 


6.06+0.14(4) 


7.20+0.33 (9) 


7.18+0.21 (5) 


7.11+0.18(19) 


35468 


y Ori 


0.088+0.024 (6) 


7.16: 


t0.26(l) 


5.98+0.08 (3) 


7.00+0.19(7) 


7.06+0.21 (3) 


7.23+0.14(11) 


36591 




0.059+0.016 (9) 


7.33: 


t0.20 (1) 


6.21+0.15 (3) 


7.25+0.29 (10) 


7.06+0.21 (1) 


7.32+0.19(31) 


50707 


15 CMa 


0.100+0.024 (7) 


7.34: 


t0.20 (1) 


6.16+0.17 (3) 


7.31+0.30 (8) 


7.17+0.21 (1) 


7.23+0.19(16) 


51309 


I CMa 


0.047+0.020 (7) 


7.30: 


t0.39(l) 


5.92+0.21 (3) 


7.18+0.31 (9) 


6.99+0.19(8) 


7.18+0.30(10) 


52089 


eCMa 


0.078+0.018 (5) 


7.39: 


1:0.23 (1) 


6.08+0.13 (3) 


7.23+0.24 (10) 


7.16+0.25 (4) 


7.16+0.15 (22) 


111 123 


/3Cru 


0.067+0.020 (9) 


7.30: 


1:0.20 (1) 


6.18+0.20 (4) 


7.25+0.23 (9) 


7.14+0.24 (1) 


7.23+0.24 (17) 



onance line. This feature is strongly blended with metallic lines 
(e.g. Mnm /12065.9) and is usually weak in the objects suffer- 
ing significant depletion. The bulk of the boron data is derived 
from Proffitt & Quigley (I2001I I who used archival International 
Ultraviolet Explorer (lUE) high-dispersion SWP spectra. They 
mention that excellent fits were achieved for T^ff < 28 000 K, 
which is the case for all stars with boron data in our sample. 
Contrary to Bii /11362, for instance, the NLTE corrections for 
B III /12066 are modest in the Teff range of interest (typically 
-0.15 dex) and should not constitute a major source of error. 
Although it is difficult to assess the reliability of these abundance 
determinations, subsequent studies based on higher quality data 
obtained with the STIS spectrograph onboard the Hubble Space 
Telescope (HST) support these estimates (Brooks et al. 120021 
Venn et al. 2002; Mendel et al. 2006). Namely, the same results 
within the uncertainties were found in cases of detection (e.g. 12 
Lac), while more stringent upper limits were set for stars with 
very weak lines (e.g. 15 CMa and /3 Cep). These more recent 
values were preferred and are those quoted in TableU] 



5. Results 

Figure [3] shows the behaviour of the [N/C] and [N/O] ratios, as a 
function of the boron content. The targets are plotted using dif- 
ferent symbols according to the detection, or lack thereof, of a 
magnetic field. Five stars lacking boron data, but which have 
been searched for magnetic fields are plotted off scale to the 
right-hand side of this figure. Three groups of stars with distinct 
chemical properties, and which may define an evolutionary se- 
quence, can be identified (see boxes in the upper, left-hand panel 
of FigO and are discussed in turn below (see also Table |4]for a 
classification of the targets in the different groups). 

Group I - First, 5 stars exhibit [N/C] and [N/O] ratios close 
to -0.6 and -0.9, respectively. These values are nearly identi- 
cal to those found for the Sun, either using classical 1-D model 
atmospheres (Grevesse & Sauval 1998) or time-dependent, 3-D 
hydrodynamical models incorporating updated atomic data and 
NLTE corrections (Asplund et al. 2006). They can thus be con- 
sidered as the baseline values in the absence of deep mixing. 
This is supported by the lack, or relatively mild, boron depletion 
compared with the meteoritic value (log e[B]=2.78 dex; Zhai & 
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Fig. 1. Variations of the CNO abundances and their ratios, as a function of the effective temperature, surface gravity and pro- 
jected rotational velocities. The dashed lines indicate the typical values found for nearby B dwarfs (M06; Daflon & Cunha 20041 
Gummersbach et al. 119981 Kilian -Montenbruck et al. I1994I I. Similar [N/C] and [N/O] ratios are found for the Sun (Grevesse & 
Sauval il998. Asplund et al. I2006i l. and can thus be considered as the baseline values in the absence of rotational mixing. 



Shaw ll994] l. These stars have not yet undergone substantial deep 
mixing and boron is only beginning to be transported to deeper 
layers of the stellar envelope where it is destroyed by proton cap- 
ture. No clear case for an N-enriched star with a solar B abun- 
dance can be found (the N overabundance of v Eri and fi Cru is 
only marginal), in accordance with the theoretical expectations. 

Group II - Two stars in'^ Ori and HD 36591) present an in- 
teresting behaviour: they are very boron depleted (by a factor of 
at least 40), yet have solar nitrogen abundances. Shallow mix- 
ing has dramatically depleted the photospheric boron, but deep 
mixing is either absent or, at least, has yet to bring detectable 
amounts of CN-cycle burning products to the surface. As can 
be seen in Fig|3] the chemical properties of these two stars are 
not well-reproduced by theoretical models (Fliegner et al. 11 9961 
Heger & Langer 2000). Specifically, the surface B depletion and 
N enrichment do not seem to proceed gradually and in parallel, 
as predicted, but that the latter is delayed and may only begin 
after boron has been depleted by ~ 1.5 dex. 

Group III - Finally, another population (8 stars) exhibits 
very low boron abundances (indeed mostly upper limits) coupled 
with an N enrichment reaching a factor about 3-4 (the separa- 
tion between the two subsamples of N-normal and N-rich stars 
is clearer in terms of the [N/C] ratio, which is the most robust 
diagnostic for an N excess). The abundance patterns of this sub- 
group are well-reproduced by the theoretical models of Heger & 
Langer (|2000 ). but only for initial velocities exceeding -200 km 



s"' (see Fig|3j3 This is well above the true rotation rates of (, 
Cas, 5 Get and (5 Cep: il/? < 60 km s ^ The other N-rich stars 
have ilRsin/ < 50 km s ' (Table |3]l. We will argue below that 
these stars only suffered a moderate amount of angular momen- 
tum loss along the main sequence and that they were likely to be 
akeady slow rotators on the ZAMS. 



6. Discussion 

6.1. A population of slowly -rotating, N-rich Galactic B dwarfs 

Our NLTE abundance analysis confirms the results of previous 
studies (e.g. Kilian 1992 ; Gies & Lambert 1 1992| l and supports 
the existence of a population of B stars which, despite having 
low rotation rates, displays the observational signature of deep 
mixing during the early phases of their evolution. Of the 36 
(non-supergiant) B stars observed by Gies & Lambert ( |1992| l, 
about one third exhibit a [N/C] ratio above about -0.4 dex. Such 
a high incidence of N enriched stars in a primarily magnitude 
limited (with an additional constraint on the projected rotational 
velocity) sample suggests that such objects make up a significant 
fraction of all B dwarfs in the solar neighbourhood. A sizeable. 



^ The Geneva models do not give predictions for the boron abun- 
dances, but they lead to even lower N enrichments (Meynet & Maeder 
|2003t . 
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Fig. 3. [N/C] and [N/O] ratios as a function of the NLTE B abundances: (a) y Peg, (b) f Cas, (c) 6 Cet, (<i) v Eri, (e) tt'* Ori, (/) y 
Ori, ig) HD 36591, (/i) /3 CMa, (/) CMa, (j) 15 CMa, (A;) t CMa, (0 e CMa, (m) HD 85953, (n) y6 Cru, (o) V836 Cen, (p) r Sco, 
(q) 9 Oph A, (r) V2052 Oph, (i) yS Cep and (f) 12 Lac. Filled circles: stars with a detected magnetic field; open circles: stars with 
null detections; crosses: stars without magnetic data (in the online version of this journal, red and blue circles refer to stars with and 
without a detected magnetic field, respectively). Five stars lacking boron data, but which have been the subject of magnetic searches 
are shown off scale to the right-hand side of this figure. The arrows denote upper limits. The abundance data are compared with the 
theoretical predictions of Heger & Langer ^20001) for a 12 M© star and three different values of the equatorial rotational velocity 
on the ZAMS: 0/?zams=99 {left-hand panels), 206 {middle panels) and 328 km s"' {right-hand panels). The locus in each panel 
{dashed line and open squares) defines an age sequence with the time elapsed from the ZAMS increasing leftwards: from f=0 to 
15 Myrs (0 to 12.5 Myrs for QRzkm5-99 km s"') in steps of 2.5 Myrs (see upper, right-hand panel). The initial [N /C], [N /O], and 
boron abundances at t-Q have been taken as -0.6, -0.9 (representative of the pristine solar value: Grevesse & Sauval [T998i Asplund 
et al. 120061) and log e(B)=2. 78+0.04 dex (meteoritic value; Zhai & Shaw 1994), respectively. The boxes in the upper, left-hand panel 
delineate the three classes of stars (Groups I, II and 111) with different chemical properties (see text). 



similar population has also been shown to exist in the Magellanic 
Clouds (Hunter et al. l2007al b; Trundle et al. l2007b . 

Mass transfer processes in close interacting binaries may 
also dramatically alter the CNO surface abundances of the two 
components (e.g. Vanbeveren il989j and may be responsible for 
the existence of the so-called OBN stars whose unusual chemical 
properties are believed to be connected to binarity (e.g. Levato 
et al. 119881 ). However, one expects in that case a much higher 
N overabundance and C depletion than observed in our sample 
(Schonberner et al. 11988 1). A good example of such a star with 
a more extreme CNO abundance pattern is the blue straggler 6 
Car (Hubrig et al., in prep.). On the other hand, fi Cep is the only 
N-rich star known to be member of a binary system. The orbit is 
very wide ('Porb~90 yr; Pigulski & Boratyn [19921 1 and the exis- 
tence of a past episode of mass transfer is debatable. In the case 
of /r"* Ori, the strong boron depletion but roughly normal nitro- 
gen abundance have been claimed to be the unambiguous sig- 
nature of rotational mixing (Fliegner et al. 1996). Furthermore, 
three other single-lined binary systems (y Peg, fi Cru and 9 Oph 
A) have normal abundances. Summarizing, mass exchange pro- 
cesses do not appear as a viable explanation for all the N over- 
abundances observed. 

With the caveat that not all the N-rich stars in our sample 
have been the subject of intensive monitoring campaigns, one 
can note that at least two stars in this subsample (y Ori and e 
CMa) are classified as poor or rejected y6 Cephei candidates in 
the recent catalogue of Stankov & Handler (120051 ). As can be 
seen in FigHl two N-rich stars (( CMa and r Sco) also fall well 



outside the instability domains for pressure and gravity modesQ 
On the other hand, some well-known pulsators are N-normal 
(e.g. y Peg, (3 CMa). This argues against a purely pulsational 
origin for the nitrogen enrichment. 

On the other hand, except perhaps in the case of i CMa, 
the stellar winds are too weak to remove the outermost stellar 
layers and expose boron-poor material. As an illustration, mass- 
loss rates M~10"^ M© yr"' are required (e.g. Venn et al. 2002]l, 
whereas two of the most evolved stars in our sample fall short by 
at least one order of magnitude: M~6 and 9 x 10"'^ Mq yr"' for 
P CMa a nd e CMa, respectively (Cassinelli et al. |1995l Drew et 
al. 119941 ). The same conclusion holds for the main-sequence star 
T Sco, with M~2 X 10"** Mq yr"' (Donati et al. 2006b). 

It can be argued that some stars in our sample were rapid ro- 
tators on the ZAMS, hence experiencing strong rotational mix- 
ing and surface N enrichment, but that substantial loss of angu- 
lar momentum ensued, spinning down the stars to the observed 
levels. This hypothesis is not supported by evolutionary mod- 
els, which only predict a slow decline of the rotation rate during 
core-hydrogen burning (e.g. Meynet & Maeder 120031) . as con- 
firmed by observations of early-B dwarfs in young open clus- 
ters (e.g. Huang & Gies 2006; Wolff et aX.IMT). However, the 
case of the magnetic stars is more controversial. For instance, 

^ A blue loop evolution can be safely ruled out for the evolved star t 
CMa according to evolutionary models (both with and without rotation 
included). The observed abundance anomalies are indeed much less se- 
vere than expected if the star had gone through the first dredge-up (e.g. 
Venn [T999] >. 
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Fig. 4. Position of the programme stars in the Hertzsprung-Russell diagram: (a) y Peg, (h) ^ Cas, (c) 6 Cet, (d) v Eri, (e) n'^ Ori, 
(f) y Ori, (g) HD 36591, (h) p CMa, (/) CMa, (j) 15 CMa, ik) l CMa, (/) e CMa, (m) HD 85953, in) /3 Cru, (o) V836 Cen, (p) r 
Sco, iq) 6 Oph A, (r) V2052 Oph, (s) (3 Cep and (f) 12 Lac. 0/?en circles: Group I (B- and N-normal stars), open squares: Group 
II (N-normal, but B-depleted stars), filled circles: Group III (B-depleted and N-rich stars). See Fig[3]for a definition of the three 
subgroups. In the online version of this journal, green and magenta symbols refer to stars with and without an N excess, respectively. 
The stars lacking boron data (see right-hand side of FigO are included in this figure assuming that the 3 N-normal and 2 N-rich 
stars belong to Groups I and III, respectively. The size of the symbols is linearly proportional to the OR sin / values (see Table [3]l. 
The luminosities were computed using Hipparcos parallaxes and the bolometric corrections of Flower (il996i) . The extinction in 
the V band. Ay, was derived from the theoretical {B - V) colour indices of Bessell et al. ( 119981 ). The evolutionary tracks for solar 
metallicity and without rotation are taken from Schaller et al. (1992). The initial masses are indicated on the right-hand side of 
this figure and the ZAMS is shown as a dotted line. The theoretical instability strips for fi Cephei and SPB stars are shown as long 
dashed and dotted-dashed lines, respectively (Miglio et al. 120071 ). The models have been computed assuming a metallicity Z=0.01 
(value typical of the fi Cephei stars in our sample; see M06), updated OP opacities, the solar mixture of Asplund et al. (i2006j and 
the mean Ne abundance of Cunha et al. ( I2006I I. 



a magnetically-confined wind could constitute an efifective way 
of dissipating angular momentum (as in the case of the O star 
HD 191612; Donati et al. |2006al ). The influence of a large-scale 
magnetic field with a dipole configuration on a radiatively-driven 
outflow can be parametrized by a dimensionless confinement pa- 
rameter. 



(3) 



where Bpoi is the polar field strength, R-i, is the stellar radius, M 
is the mass-loss rate and Vco is the terminal wind velocity (ud- 
Doula & Owocki'2002). The emergent flow is channeled along 
closed magnetic loops and deflected towards the magnetic equa- 
torial plane for rj > I. Assuming an oblique rotator model for 
^ Cas, V2052 Oph and /3 Cep leads to a polar field of about 
300 G (Neiner et al. |2003al b; Henrichs et al. |2007l ). The field 
geometry is more complex in r Sco, but the average field emerg- 
ing at the surface typically reaches similar values (Donati et al. 
|2006bl . If we adopt, as representative values, Bpoi=300 G, R*-5 
Ro, M=10"** M© yr^' (see Sect. [63 and v„ = 1500 km s^', we 



obtain 77~30. It is thus likely that the stellar outflow would be 
magnetically confined within the Alfven radius, Ra, in our mag- 
netic targets. Following Donati et al. ( 1200 11 1, the magnetic brak- 
ing timescale, r, can be roughly estimated through the following 
expression 



(R^ 



m \R 



(4) 



where k is the fractional gyration radius, is the stellar mass 
and m is the effective mass-loss rate (i.e. the amount of wind 
material not trapped within the magnetosphere). Characteristic 
values for m and Ra have been obtained in the case of /? Cep 
and T Sco by Donati et al. (2001 2006b) using the magnetically- 
confined wind-shock model of Babel & Montmerle ( 1997 ). This 
led to magnetic braking timescales that are much longer than the 
estimated stellar ages (1 10 Myrs and 5 Gyrs, respectively). There 
is considerable uncertainty surrounding m and the loss of angu- 
lar momentum through magnetic phenomena is not fully under- 
stood (for instance, some chemically peculiar Bp stars undergo 
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Fig. 2. Variations of the He, Mg, Al, Si, S and Fe abundances, 
as a function of the effective temperature. The mean values are 
indicated by dashed lines. 



unexpectedly strong, sometimes abrupt, increases in their rota- 
tional period: Mikulasek et al. 2007 ; Trigilio et al. |2007| l. but 
it appears unlikely that the slow rotation rate of the N-rich stars 
results from magnetic braking. Another argument against this in- 
terpretation is the fact that some non-magnetic stars in our sam- 
ple are also intrinsically slow rotators (Table O. This suggests 
that another mechanism could also be at work, a possibility be- 
ing that this property is inherited from the pre-main sequence 
phase (e.g. St^pieril^OH Wolff et al. .2007j . 

6.2. A link with magnetic fields? 

Rotational mixing effects are expected to be more conspicuous 
at lower metallicities and for increasing stellar mass, age and ro- 
tation rate (e.g. Meynet & Maeder [20001 ). There are no trends 
between the abundance data and OR sin / (Fig[T]i or O/?, but this 
is somewhat expected considering the limited range in rotation 
rate sampled by our sample (and the noise introduced by the un- 
known inclination angles in the case of the projected rotational 
velocities). More telling is the lack of evidence for the stars 
with an N excess to be on average older and/or more massive 
than the stars with normal N abundances (FigHji. As can also 



be seen in Fig|3] all the N-rich stars reach very similar ("sat- 
uration"?) levels of N enhancement (~0.5 dex) despite having 
different physical properties. The inadequacy of the theoretical 
models to reproduce our observations is well-illustrated by the 
unexpected detection of a factor of 3 nitrogen excess in r Sco 
and V2052 Oph. First, these stars are very slow rotators (espe- 
cially T Sco with QJi~6 km s Donati et al. 2006b) and rota- 
tional mixing should be particularly inefficient in that case. This 
problem could be lessened if one assumes a much steeper inter- 
nal rotation law than cuiTently assumed, leading to an increased 
transport of the chemical elements through shear instabilities. A 
decline of the rotational velocity by a factor 3-5 from the core 
to the surface has been inferred from asteroseismological stud- 
ies in V Eri (Pamyatnykh et al. I2004I I and V836 Cen (Dupret et 
al. l2004l l. Although this information is unfortunately lacking for 
the N-rich stars, this appears in accordance with the predictions 
of evolutionary models (Meynet & Maeder 2000). Second, these 
stars are also probably very young (see FiglU and should not 
have had time to build up a significant N overabundance at the 
surface. Taken together, this suggests that an additional physical 
process to rotational mixing might play a key role in the appear- 
ance of N-enriched material at the stellar surface. 

In this respect, it is remarkable that all observational data 
collected to date point to a higher incidence of a magnetic field 
in stars with peculiar abundances: all observations of the N-rich 
stars led to a magnetic field detection, whereas all stars with 
solar abundances remained undetected^ The only two excep- 
tions to this rule are y Peg and HD 85953 (Fig|3]l. The /3 Cephei 
star y Peg has been claimed to possess a very weak magnetic 
field with a longitudinal component varying from -10 to 30 G 
(Butkovskaya & Plachinda ,2007 ). On the other hand, a magnetic 
field has been detected in the Slowly Pulsating B (SPB) star HD 
85953 at the 3-cr level in 3 out of 5 measurements spread over 
3.5 yrs with FORSl/VLT (Hubrig et al. 2006; Hubrig et al., in 
prep.). The latter example illustrates the paucity of the magnetic 
fields measurements in B-type stars and that caution should be 
exercised when using these data. The field detections reported 
in Table [T] are for the most part at the limit of the instrumental 
capabilities (especially with the older generation of spectropo- 
larimeters such as MuSiCoS) and are based on snapshot obser- 
vations with a poor temporal sampling. Only in the cases of ^ 
Cas, T Sco, V2052 Oph and /3 Cep have sufficient time-resolved 
measurements been collected to provide detailed information on 
the field strength and global morphology (Neiner et al. 2003a b; 
Donati et al. I2006bl Henrichs et al. l2007l l. Considering that the 
null detections are also often based on a handful of measure- 
ments and might be refuted when more intensive and precise ob- 
servations are undertaken, the separation between magnetic and 
non-magnetic stars should only be regarded as arbitrary at this 
stage. It is also conceivable that the N enrichment took place 
when the star was very young and going through a magnetic 
phase (similarly to 0' Ori C), and that the field strongly de- 
cayed afterwards (e.g. Hubrig et al. 2007) to a point where it 
is no longer detectable with current instrumental facilities. In 
this context, it is interesting to note that an N excess is already 
present in the young magnetic stars t Sco and V2052 Oph (see 
FiglDi. Although highly speculative, this might explain why a 



We do not discuss here the case of the N-rich, apparently slowly- 
rotating yS Cephei star HD 180642 (Morel & Aerts 2007), as the mag- 
netic field detection is not completely secure and needs to be confirmed: 
although several Zeeman features are seen in Stokes V spectra obtained 
with FORSIA'LT, none of the two spectropolarimetric measurements 
formally exceeds the Scr detection threshold (Hubrig et al., in prep.). 
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slowly-rotating, N-ric h dwa rf such as 1 Cas ([N/C]=-0.07±0.20 
dex; Gies & Lambert '1992') remained undetected during the in- 
tensive MuSiCoS observations of Schnerr ( 120071 ). A correlation 
between the amount of CN-cycle processed material at the sur- 
face and the field strength would strongly support a connec- 
tion between deep mixing and magnetic phenomena. There is 
no trend in our data between [N/C] and the rms longitudinal 
field (Bi) (Table [TJ, but it should be kept in mind that the ro- 
tational period is poorly sampled in most cases and that the 
magnetic observations are restricted to longitudinal field mea- 
surements, which are aspect dependent, and to fields larger than 
~40-50 Gauss. Despite these limitations, however, forthcoming 
magnetic observations of OB stars conducted using the new gen- 
eration of spectropolarimetric instruments (such as FORS 1/VLT, 
ESPaDOnS/CFHT or Narval/TBL) will clearly prove essential 
to establish the existence of a link between a nitrogen enrich- 
ment and magnetic fields. Such a relationship could in principle 
help to select the most promising targets for magnetic field sur- 
veys of massive stars (indeed, an N excess was already one of 
the selection criteria used by Schnerr 2007}. 

Strong, ordered kG-strength magnetic fields have been de- 
tected in two slowly-rotating 0-type stars, namely in the 
06.5f?pe-08fp star HD 191612 (Donati et al. |2006a') and in 
the young, peculiar 07 star 0' Ori C (Wade et al.[2t)06). To our 
knowledge, however, no quantitative CNO abundance analysis 
has been performed for either of these two objects. Only in the 
former has a nitrogen enhancement been suggested based on the 
observed line strengths (Walborn et al. 120031 ). Unfortunately, a 
similar lack of CNO abundance data also applies to classical, 
magnetic Bp stars lying in a similar temperature domain as our 
programme stars. Two strongly magnetic. He-rich dwarfs (HD 
66522 and HD 96446) have been analyzed under the assump- 
tion of LTE by Zboril & North ( 1999). The CNO abundances do 
not clearly differ from those of the non-magnetic stars in their 
sample, but no firm conclusion can be drawn from such a lim- 
ited number of objects. Systematic studies have been conducted 
for Ap stars, both for boron (e.g. Leckrone l 19801 1 and CNO (e.g. 
Roby & Lambert 1990). However, these stars are much cooler 
than our targets and diffusion effects have long been known to 
play a key role in the dramatic surface abundance anomalies ob- 
served in this temperature range. 

Nevertheless, radiatively-driven microscopic diffusion has 
also been proposed (Bourge et al. 120071 ) as a possible explana- 
tion for the appearance of core-processed material at the surface 
of slowly-rotating /3 Cephei variables (we recall that 5 N-rich tar- 
gets are confirmed members of this class; see Table [3]). Support 
for the occurrence of diffusion comes from the fact that it may 
also account for the accumulation of iron in the driving zone 
and the excitation of the pulsation modes through the k mecha- 
nism (e.g. Bourge et al. 120061 ). Magnetic fields were not taken 
into account in these models, but they are expected to amplify 
the transport of the ionized species along the field lines and to 
lead to a large-scale, inhomogeneous distribution of some ele- 
ments at the stellar surface. However, the Mg, Al, Si, S and Fe 
abundances do not reveal any clear deviation from a scaled solar 
pattern (see Fig|2]) and do not provide evidence for the existence 
of diffusion effects capable of significantly altering the surface 
abundances. This is supported by a comparison of the abundance 
properties of the two subsamples of N-normal and N-rich stars, 
which reveals no significant differences for chemical elements 
other than nitrogen (Table |6]). Two stars apparently have a low 
helium content (HD 36591 and t CMa), but this is at odds with 
the expected/observed position of the He weak stars in the T^ff- 
log g plane (Babel il996i Bychkov et al. 120031 ) and may be re- 



Table 6. Comparison between the chemical properties of the two 
subsamples of N-normal and N-rich stars. 





N-normal sample 


N-rich sample 




(10 stars) 


(10 stars) 


He/H 


0.071 ±0.007 


0.087±0.020 


loge(C) 


8.21±0.09 


8.11+0.10 


loge(N) 


7.67±0.11 


7.98+0.09 


loge(O) 


8.51±0.09 


8.44+0.16 


log 6(Mg) 


7.39±0.11 


7.38+0.12 


l0g6(Al) 


6.13±0.08 


6.10+0.12 


l0g6(Si) 


7.18±0.07 


7.19+0.10 


loge(S) 


7.20±0.11 


7.15+0.11 


log 6(Fe) 


7.28±0.10 


7.26+0.08 



garded with some suspicion. Slightly subsolar He abundances 
are in general derived within our sample (Fig|2]l. We suspect that 
the adopted microturbulent velocities derived from the analysis 
of the 11 features are not adequate for properly modelling the 
He I lines, and that lower values should ideally be employed in- 
stead (as also concluded for Nii; Sect.|4|i. 

Recent codes attempting to model rotational mixing in the 
presence of magnetic fields (Maeder & Mevnet [20051 Heger et 
al. 2005) also face difficulties for accounting for the data pre- 
sented in this paper. For instance, they predict a noticeable he- 
lium enrichment, which is not observed (the only exception is 
V2052 Oph, but the slight He excess is more likely to be due 
to magnetic phenomena; see Neiner et al. '2003b and M06). 
Furthermore, these models assume a magnetic field created via 
dynamo action in the stellar envelope (Spruit I2002I I. whereas 
at least three N-rich stars in our sample (^ Cas, "^^^2052 Oph 
and /3 Cep) have been claimed to possess a global dipole field 
most likely of fossil origin (Neiner et al. 2003a,b; Henrichs et 
al. 2007) . A remnant field from the star formation phase is also 
more likely in t Sco, although the field topology is probably 
much more complex (Donati et al. '2006b). In all cases, the peri- 
ods found in the magnetic data are compatible with those mod- 
ulating the UV wind line profiles, thus supporting an oblique 
rotator model. Clearly much progress remains to be done before 
the production of magnetic fields in massive stars and the mixing 
processes taking place in their interiors are satisfactorily under- 
stood (see, e.g. Zahn et al. 120071 ). 
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